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Abstract 

We report on an exclusive search for the X(3872) state in the decay 7r + 7r"J/^, J/tp — ► l + U (I = 
e,/x) from untagged 77 fusion and initial state radiation production using a 15.1 fb _1 data sample 
with the CLEO III detector. By taking advantage of the unique 7r + 7r~J/V> decay kinematics, 
separate measurements for each production process are obtained. No signals are observed and 
preliminary upper limits have been established as (2J + l)r 77 (X(3872)),B(X(3872) — > ir + ir~ J/VO 
< 12.9 eV and r cc (X(3872))fi(X(3872) -► tt+tt-J/^) < 8.0 eV, both at a 90% confidence level. 



*Submitted to the 32 nd International Conference on High Energy Physics, Aug 2004, Beijing 
T On leave of absence from University of Chicago. 
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The Belle Collaboration reported the observation of a narrow state, X(3872), in the 
decay of 152 million BB events through B^ — > K ± X, X — > n + ix~ J/ip, J/ift — > l + l~ (I = e,/i) 
measuring M(X) = 3872.0 ± 0.6 (stat) ± 0.5 (syst) MeV/c 2 and T < 2.3 MeV/c 2 at a 90% 
confidence level (C.L.) The observation was confirmed by the CDF II Collaboration 
and D0 Collaboration in inclusive production with X(3872) decaying to 7r + 7r - JM, 
J/ip — > fi + fi~ in pp collisions at y/s = 1.96 TeV and also by the BaBar Collaboration |4( 
using 117 million BB events with the same decay channels as the Belle observation. CDF 
II, D0, and BaBar measured masses of 3871.3 ± 0.7 (stat) ± 0.4 (syst) MeV/c 2 , 3871.8 ± 
3.1 (stat) ± 3.0 (syst) MeV/c 2 ', and 3873.4 ± 1.4 (stat) MeV/c 2 , respectively, and widths 
consistent with their respective detector resolutions. 

A large number of theoretical papers now exist BQBHHHQQQElIElI, 



171 . Il8j with different interpretations of the nature of the X(3872) state and its possible 



quantum numbers. Among the possibilities are: (a) a charmonium state; (b) a loosely 
bound "molecular" state; or (c) an exotic state. Barnes and Godfrey ffj and Eichten, Lane, 
and Quigg 0| have examined the charmonium options in detail and conclude that, because 
of the small width of X(3872) and despite the larger predicted masses of the candidate 
charmonium states, the viable charmonium options are the negative C parity states of 1 3 D2 
(jpc = 2~), 1 3 D 3 (3~), and 2 1 P 1 (1+-), and positive C states of 1 X D 2 (2~+) and 2 3 Pi 
(1 ++ ). Tornquist and Swanson [l(| propose that, since X(3872) is close to the DD* 
threshold (M(D°) + M(D *) = 3871.3 ± 0.5 MeV/c 2 [jIj) and there are no published studies 
in 7r°7r°J/?/> decays, the observed decay 7r + 7r~J/?/> proceeds through p°J/ip, p° — > n + 7i~ . They 
suggest h and 1 ++ as the most likely assignment in a molecular model. Finally, we note 
that Morningstar and Peardon [Tjj predict a three gluon 1 glueball with a mass of 3850 
MeV/c 2 . 

Searches for X(3872) have been made but no positive signals have been observed in 
the decay channels X(3872) -> 7 x c i, 7Xc2, rjJ/ip, D+D", D°D°, and D°D 7r° 0, H 

2l|,|23. C. Z. Yuan, X. H. Mo, and P. Wang determined an upper limit for initial state 



radiation production of r ee (X(3872))£(X(3872) -> tt+tt-J/^) < 10 eV (90% C.L.) |23j using 
a previously published 22.3 pb -1 data sample of e + e~ annihilations at y/s = 4.03 GeV with 
the BES detector j2J|. With all of these searches, the new state X(3872) has only been 
experimentally observed decaying to 7r + 7r _ J/^. 

The variety of possibilities for the structure of X(3872) suggests that, irrespective of the 
models, it is useful to limit the J PC of X(3872) as much as possible. Untagged 77 fusion, 
which only populates states of positive C parity and even values of total angular momentum, 
can be used to shed light on the charge parity and total angular momentum of X(3872) and 
its population via initial state radiation can determine its leptonic width and its likely vector 
character. The present investigation is designed to study the possible production of X(3872) 
in untagged 77 fusion and initial state radiation in the final state 7r + 7r~J/?/>, J/ip —>■ (/ 
= e,/i). 

Untagged 77 fusion and initial state radiation (ISR) resonance production have similar 
event characteristics. Untagged 77 fusion resonance production occurs when a photon is 
emitted by each incident electron/positron, the two almost real photons "fuse" to form the 
resonance, while the scattered electron/positron are not detected. ISR resonance production 
predominately occurs when either the initial electron or positron emits a hard photon which 
lowers the center-of-mass energy of the electron-positron system to the invariant mass of 
the produced resonance before annihilation. The radiated photon also has an angular dis- 
tribution which is very sharply peaked along the beam axis and therefore is rarely detected. 
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Both types of resonance production processes have a total observed energy (E tota i) much 
smaller than the center-of-mass energy of the electron-positron system and a small observed 
transverse momentum. 

The data used for this X(3872) search were collected with the detector in the CLEO III 
configuration [25| at the Cornell Electron Storage Ring (CESR). The CLEO III detector is a 
cylindrically symmetric detector designed to study e + e~ annihilations and provide 93% cov- 
erage of solid angle for charged and neutral particle identification. The detector components 
important for this analysis are the drift chamber (DR), Csl crystal calorimeter (CC), and 
muon identification system (MIS). The DR and CC are operated within a 1.5 T magnetic 
field produced by a superconducting solenoid located directly outside of the CC. The DR 
detects charged particles and measures their momenta and ionization energy loss (dE/dx). 
The CC allows precision measurements of electromagnetic shower energy and position. The 
MIS consists of proportional chambers placed between layers of the magnetic field return 
iron to detect charged particles which penetrate a minimum of three nuclear interaction 
lengths. 

The data consist of a 15.1 fb -1 sample at or near the energies of the T(nS) resonances, 
where n = 1-5, and in the vicinity of the AbAb threshold. Table 1 lists the six different 
initial center-of-mass energies (CME) and integrated luminosities of the data sample. 

TABLE I: Data sample used for this X(3872) search, (i/s ) is the integrated luminosity aver- 
aged initial CME for the T(1S-5S) resonance and threshold AbAb data samples, respectively, and 
£(e + e~) is the e + e~ integrated luminosity for each CME region. 





<Vs ) (GeV) 


£(e+e-) (fb- 1 ) 


T(1S) 


9.458 


1.47 


T(2S) 


10.018 


1.84 


T(3S) 


10.356 


1.67 


T(4S) 


10.566 


8.97 


T(5S) 


10.868 


0.43 


AbAb threshold 


11.296 


0.72 



The expected characteristics for 77 fusion and ISR resonance production are studied by 
generating signal Monte Carlo (MC) samples. The 77 fusion resonance signal MC process 
is e + e~ — > e + e~(77), 77 — > X(3872) using the formalism of Budnev et al. The ISR 

signal MC process is e + e~ — > 7X(3872) using the angular distribution expression from M. 
Benayoun et al. [28j]. The X(3872) and are decayed according to phase space. The MC 
simulation of the CLEO III detector was based upon GEANT 3.21/11 [26[ with simulated 
events processed in the same manner as the data. 

A fully reconstructed event has four charged particles and zero net charge. All charged 
particles must lie within the DR volume and satisfy standard requirements for track quality 
and distance of closest approach to the interaction point (IP). Events must also have detected 
Etotai < 6 GeV. Signal-to-background studies are performed to optimize signal efficiency and 
background suppression using the 77 fusion signal MC sample for the signal and background 
from data with M(7t + tt~1 + 1~) — M(/ + Z~) above and below the X(3872) search region, i.e., 
M(7r + vr-/ + /-) -M(Z + /~) = 0.63-0.7 or 0.85-0.92 GeV/c 2 . The selection variables optimized 
are the total neutral energy (E neu ) of the event, total transverse momentum of the four 
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FIG. 1: AM = M(vr+7r-/+r) - M(Z+r) data distribution for AM = 0.514-0.850 GeV/c 2 . 
ip(2S) is clearly visible and there is no enhancement in the X(3872) region. 
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charged tracks (pt r ), lepton pair invariant mass (M(/ + Z~)) of the J/ip — » decay, and 
particle identification (PID) of the charged tracks. Based on the optimization studies, events 
are selected with E neu < 0.4 GeV and p tr < 0.3 GeV/c. Events with a J/ip — > e + e~ decay 
have both electron candidates satisfying dE/dx and shower energy criteria consistent with 
the electron hypothesis and an invariant mass M(e + e~) = 2.96-3.125 GeV/c 2 . Events with 
a J/ip — * decay have both muon candidates appearing as a minimum ionizing particle 

in the CC, at least one muon must penetrate a number of interaction lengths in the MIS 
consistent with its momentum, and an invariant mass M(/i + /i _ ) = 3.05-3.125 GeV/c 2 . Each 
of the two pions recoiling against the J ftp needs to individually satisfy the dE/dx pion 
hypothesis. 

Fig. 1 shows the AM = M(7r + 7r~/ + /~) - M(/ + /~) data distributions for events which 
pass the selection criteria and have AM = 0.514-0.850 GeV/c 2 . A ip(2S) signal is clearly 
visible while no enhancement is apparent in the X(3872) region. A MC sample of ip(2S) 
produced via ISR is generated to determine the efficiency for detecting the observed ■0(2S) 
signal, from which an expected number of ip(2S) — > 7r + 7r~ J/ip, tt + tc~(1 + 1~) events is derived 
to be 226 ± 11 events, where the error is from the MC efficiency and J/ip — ■> / + /~ PDG 
branching fraction [if uncertainties. The observed number of ip(2S) events is determined by 
fitting the ^(2S) region with the detector resolution and a mass-independent background. 
The detector resolution is determined by fitting the ip(2S) ISR MC resolution function with 
a double Gaussian. The narrow Gaussian width and relative area and width of the double 
Gaussian are fixed when the ^(2S) is fitted. The observed number of ip(2S) is 206 ± 15 
events and consistent with expectations. 

A feature unique to ISR production of charmonium-like resonances which decay via 
7T + 7T~J/ip at initial energies of this data sample is the correlation between the tracks of 
the two leptons. Fig. 2 shows the two-dimensional cos(#) distributions for lepton tracks in 
the X(3872) ISR resonance and 77 fusion signal MC samples. A parabolic cut is applied to 
the two-dimensional lepton pair cos(6 l ) distribution to separate the two phenomena and is 
shown in Fig. 2. This separation removes ~14% of the X(3872) 77 MC fusion signal events 
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FIG. 2: Two dimensional lepton pair cos(#) distribution for the X(3872) ISR (left) and 77 fusion 
(right) resonance MC samples. The lines indicate where the ISR resonance and 77 fusion samples 
were separated. 

but rejects more than 99.5% of the ISR events. 

The total efficiencies in the six different initial CME regions are listed in Table 2. The 
efficiencies are determined from the X(3872) 77 fusion and ISR resonance signal MC samples. 
The same selection criteria is applied to both MC samples except with their appropriate two- 
dimensional lepton pair cos(#) correlation. 

TABLE II: Detection efficiencies determined from the signal 77 fusion and ISR resonance MC 
samples. e ee (e^) is the total efficiency for detecting X(3872) in events with a 7r + 7r~J/^, J /if; — > 
e + e~ (J/V — > decay. The appropriate 77 fusion/ISR separation is applied to the respective 

signal MC sample. 







77 Fusion 


ISR 




(y/B ) (GeV) 


fee 






(■mi 


T(1S) 


9.458 


0.128(4) 


0.160(4) 


0.065(3) 


0.083(3) 


T(2S) 


10.018 


0.121(3) 


0.151(4) 


0.054(2) 


0.062(3) 


T(3S) 


10.356 


0.115(3) 


0.137(4) 


0.042(2) 


0.043(2) 


T(4S) 


10.566 


0.123(4) 


0.145(4) 


0.0186(14) 


0.0165(13) 


T(5S) 


10.868 


0.113(3) 


0.139(4) 


0.0025(5) 





AbAb threshold 


11.296 


0.104(3) 


0.126(4) 


0.0001(1) 






The AM data distributions in the X(3872) search region for 77 fusion and ISR resonance 
production are shown in Fig. 3. The number of observed events (N Q b s ) is determined by 
a maximum likelihood fit of the AM data distribution using a production-mode-dependent 
detector resolution and mass-independent background. The respective detector resolutions 
are determined by a double Gaussian fit of the signal MC resolution functions and their 
shapes are illustrated in Fig. 3. For each resolution function, the mean is fixed to the 
Belle observation of AM = 0.7751 GeV/c 2 with the narrow Gaussian width and the double 
Gaussian relative area and width also fixed. The upper limits on the observed number of 
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FIG. 3: 77 fusion (top) and ISR (bottom) AM data distributions with AM = 0.7-0.85 GeV/c 2 
satisfying respective selection criteria. The corresponding MC-determined resolution functions are 
illustrated with the means fixed at AM = 0.7751 GeV/c 2 . 



X(3872) events in untagged 77 fusion and ISR resonance production are N obs < 2.36 and 
N obs < 2.1 at 90% C.L., respectively. 

The X(3872) production results need to incorporate the six different CME regions of 
the data sample. The cross section for a 77 fusion or ISR produced resonance of mass m R 
decaying to 7r + 7r~ J /if), J /if) — > l + l~ is 

a(R) 77 , ISR £(R - 7T + 7T- J/lf)) = — — — — — (1) 

where f« = £j(e + e~)a(sj,m R ) 77sI g R is the flux for 77 fusion or ISR production, £j(e + e~) is 
the e + e~ integrated luminosity, and e ee> j (e Wi j) is the detection efficiency for events with 
7r + 7r~ J/-0, J/-0 — ► e + e~ (J/-0 — > decays at an initial CME, ^/sT. £>(R — > 7r + 7r~J/'?/>) is 

the R — » 7r + 7r~ J/-0 branching fraction while 13 (J /if) — ► e + e~) and B(J/if) — > ^ + Ml) are the 
PDG branching fractions |5[. For a narrow resonance, <r(sj,m R ) 77 is defined as [27] 

*(*.**)<n = (-£-) 2 K^)(M^) - i) 2 - s(M Jr)) 8 ]. (2) 
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and a(si,m R )isR is @ 



a(s,,m R ) ISR = J2L(21n^ - 1)(1 -x, + f ) (3) 
7rXjSj m e 2 

2E 

where Xj is the ratio of the ISR photon energy to the initial beam energy (—^J and m e is 
the electron mass. <t(R) 77 is related to the 77 fusion resonance production by 

= M2J + i)r 7T (R) 

m R 

and ct(R)isr is related to the ISR resonance production by 

12vr 2 r cc (R) 

0\R)isr = (5) 

hir 

where T 77 (R) (r ee (R)) is the two-photon (e + e~) partial width of the resonance. 

The sources of systematic uncertainty arise from possible biases in the detection efficiency 
and estimated background level. These are studied by varying the track quality, IP, 77 
fusion/ISR separation, and selection criterion optimized in the signal-to-background studies. 
Other systematic uncertainties are from the e + e~ luminosity measurement and PDG J/ip —>■ 
l + l~ branching fractions. The total systematic uncertainty contributing to 77 fusion (ISR) 
resonance production is 18.5% (23.4%), determined by summing individual contributions 
in quadrature, and is incorporated into the final results by increasing the measured upper 
limits by the respective systematic uncertainties. The specific X(3872) angular distributions 
for a given J PC have not been considered in this analysis. 

Untagged 77 fusion resonance production tests the possibility of a state having positive 
C parity and an even value of total angular momentum. The preliminary upper limit for 
X(3872) 77 fusion production is 

(2J + l)r 77 (X(3872))£(X(3872) -> J/-0) < 12.9 eV 

at a 90% C.L. Assuming i3(B ± -> K ± X(3872)) « ^(B ± -> K±^(2S)) = (6.8±0.4)xl0" 4 Q 
leads to fi(X(3872) -> tt+tt-J/^) « 0.02 for both the Belle Q and BaBar @ results. This 
translates into (2J+l)r 77 (X(38 72)) < 0.645 keV (90% C.L./ for this analysis. The X(3872) 
(2J+l)r 77 (R) is more than four times smaller than for r] c , Xco, and Xc2, i-e., (2J+l)r 77 (r7 c ) 
= 7.4 ± 2.3 keV, (2J+l)r 77 (x c0 ) = 2.6 ± 0.5 keV, and (2J+l)r 77 (x c2 ) = 2.6 ± 0.3 keV 
(l9| . respectively. If X(3872) is a vector meson, the preliminary 90% C.L. upper limit for 
ISR resonance production is 

r ee (X(3872))£(X(3872) -> tt+tt-J/V') < 8.0 eV. 

From the observed ^(2S) signal, the upper limit for the X(3872) ISR production yield is 
N X (3872)/N^(2S) < 0.011 (90% C.L.). The ISR production upper limit from this X(3872) 
search includes our estimated systematic uncertainty and is comparable to the ISR upper 
limit of r ee (X(3872))£(X(3872) -> tt+tt-J/^) < 10 eV (90% C.L.) using the BES data @. 
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